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Abstract. Disc non-axisymmetrc components, such as spirals and central bars, are nowa- 
days known to play an important role in shaping galactic discs. Here we use Tree-SPH N- 
body simulations to examine the elfect of these perturbers on two aspects: the occurrence 
of multiple patterns in discs and the elfects of radial migration on disc thickening. We 
find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop 
in all models. Interaction among these asymmetric features results in a large scale stellar 
migration. However, we show that, despite the strong radial mixing, discs cannot be thick- 
ened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as 
stars migrate radially outwards. We also find that the bulge contribution to a thick-disc 
component for an Sa-type galaxy at ~ 2.5 disc scale-lengths is less than 1% and zero in 
the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick 
disc formation via stellar radial migration. 



1 Introduction 

Major progress has been made over the last few years in understanding disc formation in cosmological 
simulations, although many discrepancies remain between the predictions of such simulations and the 
properties of the disc galaxies we observe today (e.g., ||22| ). Furthermore, there is now increasing evi- 
dence that mergers may not necessarily be the dominant process in the formation of discs, but internal 
evolution processes also play a major role in shaping galactic discs (e.g., 121151 . Recent theoretical 
studies have found that: (i) extended discs can be explained by satellite- (| 18 1) or spiral-driven migra- 
tion (|21 1); (ii) angular momentum bar-halo coupling ([JJ) and a powerful radial migration mechanism 
involving bar-spiral arms interaction ( 1141151 ') are at play in barred spirals, such as the Milky Way; and 
(iii) thick discs, traditionally explained by heating due to mergers and/or accretion of satellites ([TTl), 
could also be caused by secular evolution linked to gas-rich, turbulent, clumpy discs (|2|) or radial 
migration (|23 9|). 

It is, however, crucial to have a clear understanding of how much can be expected from secular disc 
evolution. In particular, here we would like to examine in detail the effect of radial migration on the 
formation of thick discs. To this end, we study three isolated giant galaxies from the GalMer databas43 
(1.6.51): the gSb (20% gas, u, = 220 km/s), gSa (10% gas, y,. = 280 km/s), and gSO (Identical to gSa 
but no gas) models. In these Tree-SPH N-body simulations the gSb model develops a bar similar to the 
one seen in the Milky Way (MW) ( II12I13I ). while both the gSO and gSa bars are substantially larger 
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Fig. 1. First and second columns: power spectrograms of the m = I (top) and m=2 (bottom) Fourier components 
for the gSa model. The first and second columns show time windows in the range t = [0.3 - 0.9] Gyr and 
t = [0.8 - 1.3] Gyr, respectively. Columns 3-5, top row: the differential density of the gSa disc for three different 
times, as indicated in each panel. The bars are aligned with the vertical axis and their CR and OLR radii are shown 
by the dotted red and solid blue circles, respectively. Columns 3-5, bottom row: present the incremental changes 
in angular momentum, AL, in time intervals of 200 Myr, centered on the time of the snapshots shown right above. 
Both axes are divided by the estimated rotation curve for the galaxy, thus the values displayed are approximately 
equal to galactic radius in kpc. Note that discontinuities in the spiral arms in the outer disc correspond to strong 
increases in AL. 

than that (see Fig.|2]l. We have shown (U_5|) that this has strong effects on the migration efficiency 
found in these discs, with larger bars giving rise to stronger mixing (also seen in studies of diffusion 
coefficients, e.g., M3I25I ). 



2 Multiple patterns In galactic discs 

A number of previous works have seen spiral structure in N-body or particle-mesh simulations at 
pattern speeds different from the bar's (e.g., II24I17I . see also fTj for observational evidence in the 
MW). It has been suggested that these multiple waves are coupled ( II20I10II ) and thus, longer-lived 
(> 0.5 Gvr. lfTgl ). rather than transient (on the order of a dynamical time). 

The simulations used in this work present yet another example of the propagation of multiple 
waves in galactic discs. In the first two columns of Fig. [T] we plot Fourier power spectrograms for the 
m = 1 (top, lopsided structure) and m - 2 (bottom, bar or 2-armed spiral structure) components for 
the gSa model. The time windows used for earlier (first column) and later (second column) times are 
t - [0.3 - 0.9] Gyr and t - [0.8 - 1 .3] Gyr, respectively. The bar is identifiable by the strongest feature 
in the inner disc (0-6 kpc, bottom panels); two additional features (spirals) are present in the outer 
disc at lower pattern speeds. The similarity between the Fourier spectrograms at different stages in the 
simulations (as the bar slows down) suggest that our discs also exhibits couplings between waves, with 
the bar possibly driving the spiral waves. 

To see how these multiple patterns appear in the disc morphology, in columns 3-5, top row of 
Fig. [T] we plot the face-on stellar density of the gSa model, with the azimuthally averaged component 
subtracted. As can be seen here (but also in all other models), all time snapshots show discontinu- 
ities/breaks (strong decrease in amplitude) in the spirals at the bar's outer Lindblad resonance OLR, 
as well as outside the OLR. Discontinuities in spiral arms can indicate changes in the dominant pat- 
tern or constructive/destructive interference between different spiral modes. For example, the lack of 
rotational symmetry in the snapshots in Fig.[T]is indicative of an w = 1 mode, possibly resulting from 
interference between two patterns of different multiplicity (e.g., a bar and a three-armed spiral wave). 
Strong m - 1 modes are found in the power spectrograms as well. Observationally, lopsidedness has 
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been shown to be a common property of spiral galaxies. For example, Il27l have shown that ~ 30% of 
field spiral galaxies in a magnitude-limited sample exhibit significant lopsidedness in the outer discs. 

2.1 Angular momentum exchange at spiral discontinuities 

We now study the angular momentum exchange at different radii in the gSa disc. In Fig. 2 by 1 15] we 
presented the cumulative changes in angular momentum for this same model, as a function of radius 
for the duration of the simulation (3 Gyr). In contrast, the bottom panels of rows 3-5 of Fig. [T] show 
the incremental changes in the time intervals 0.3-0.5, 0.5-0.7, and 0.7-0.9 Gyr. We estimate these as 
AL{r) - Li(r)- Lo(r), where Li(r) and Loir) are the final and "initial" (i.e., t=0.3, 0.5, and 0.7) angular 
momenta at each time step. In other words, each plot shows how much migration has occurred as a 
function of radius during the time step of 200 Myr, centered on the time of the snapshots shown in 
the first three columns of the same figure. Both axes are divided by the estimated rotation curve, thus 
the values displayed are approximately equal to galactic radius in kpc. The bar's corotation (CR) and 
OLR are shown as the dotted and solid vertical lines in each plot. Examining the changes of angular 
momentum, we note that the effect near the bar's OLR is similar to that at the CR and even stronger in 
the first time period shown. Very well defined multiple peaks in AL are present at all time steps, with 
strongest features near the spiral discontinuities. 

The relationship we find here between the changes in angular momentum and the gaps and breaks 
in the spiral structure occurring at the same radii, makes it compelling to interpret these with resonance 
overlap regions associated with the multiple patterns. Beyond the outermost break however, spiral 
features are most likely not self-gravitating because of the low densities in those regions. This is also 
suggested by the abrupt decrease in pitch angle. We nevertheless see changes at those discontinuities. 
These must be related to gravity torques |8|, rather than resonance overlap, given the long dynamical 
times at those radii. 



3 Radial migration and thick discs 

An interesting question that could help us understand the formation and evolution of disc galaxies is 
how their thick discs came into existence. Recently there has been a growing conviction that radial 
migration can result in a thick disc formation by bringing out stars with high velocity dispersion from 
the inner disc and the bulge. To test this hypothesis, here we analyze the three models above. The gSb 
model is comparable to the MW in mass, rotation curve, spiral structure strengths, and bar and bulge 
sizes. The gSa (dissipational) and gSO (dissipationless) models both develop bars and spirals stronger 
than those of the MW. Both gSb and gSO exhibit multiple patterns in their discs, similarly to gSa (see 
Sec. |2]l. As shown in Fig. [T] this gives rise to strong radial migration. To compare the three models, 
in the first column of Fig. |2] we plot the m - 2 Fourier amplitudes, A2IAQ, as a function of radius, 
where Aq is the axisymmetric component and m is the multiplicity of the pattern. These indicate the 
bi-symmetric structure in the disc. For example, for the gSO model the bar is identifiable by the smooth 
curve in the inner disc, which at later times has a maximum at ~ 3 kpc and drops almost to zero at 
~ 8 kpc. Deviations from zero seen beyond that radius are due to the spiral structure. Note that the 
spirals are strongest for the gSa model, related to the initial 10% gas fraction in its disc. The gSb model 
is the one with parameters closest to the MW: bar size ~ 3 - 4 kpc, Vc ~ 220 km/s, and a small bulge. 

To see how much discs thicken, in columns 2-3 of Fig. |2]we plot the edge-on view for each disc 
and bulge (separately) for earlier (0.25 Gyr) and later (4 Gyr for gSb and gSO, and 3 Gyr for gSa) times 
of their evolution, as indicated in each panel. Significant thickening is seen for all models, especially 
for gSa and gSO due to their larger bars. The bulges are also found to expand, not significantly so for 
the gSb. Although from these plots it is evident that the discs thicken, it is important to see how this 
is reflected in their scale-heights. These are shown in columns 4-5 at a radius ~ 2.5 disc scale-lengths, 
which is comparable to the solar distance from the Galactic center 

Surprisingly, we find that, despite the strong radial migration, the initial scale-heights, h, hardly 
double in the 3-4 Gyr of evolution for gSa and gSO, respectively, and increase by only ~ 50% in 4 Gyr 
for gSb. We look for the reason for this behavior in the next Section. 
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Fig. 2. First column: m=2 Fourier amplitudes, A2, as a function of time (color bar) for our model galaxies. 
Columns 2-3: edge-on view for each disc and bulge (separately) for earlier and later times of their evolution, 
as indicated in each panel. Columns 4-5: scale-heights, h, at ~ 2.5 disc scale-lengths. Note that in 3-4 Gyr 
of evolution with strong perturbation from bars and spirals, h only doubles despite the efficient stellar radial 
migration (see Fig.[T). 

3.1 Cooling of stellar samples migrating outwards 

Stellar samples in the inner galactic discs possess high velocity dispersions in all three components. 
It has been expected that by migrating radially outwards, stars retain their random energies and thus 
thicken the disc. If this were true, then why did we not see sufficient disc thickening even for the 
strong large-scale disc mixing seen in Fig. [T]? To answer this question, we consider the samples of 
stars used to estimate the scale-heights in Fig. |2] In the top three panels of Fig. [3]we plot the initial 
(dotted black), final (dashed blue), and net (dot-dashed red) vertical velocity dispersions, cr. ,-, cr. f, and 
Act-, respectively, as functions of the initial radius for particles ending up in the green bin at the final 
time. It is immediately apparent that stellar samples aiTiving from the inner disc decrease their velocity 
dispersion by as much as 50% when coming from the smallest radii, with the effect diminishing the 
smaller the radial distance traveled. It is also notable that the stars coming from radii exterior to the 
green bin increase their velocity dispersion. The negative slope in cr^j indicates that the outward 
migrators indeed heat the outer disc, but not significantly. As seen in the initial radial distributions at 
the bottom three panels of Fig. [3] the majority of stars populating the annuli at 2.5 disc scale-lengths, 
rd, (green bin) have arrived/migrated from the inner disc. 

Note that the fact that stellar samples migrating outwards decrease their vertical energy does not 
mean that the discs cool down. The fourth column of Fig. |3]plots the time evolution of the vertical 
velocity dispersion profiles for gSb (top), gSa (middle), and gSO (bottom). An increase in cr, is seen at 
every radius, as expected. 

To understand this cooling process we consider a simple galactic disc model, where the vertical 
action of stars, J~, is conserved as they shift guiding radii (i.e., migrate). This is true in the limit that (i) 
the vertical motion decouples from the planar motion (a good assumption when Zm-dx r which holds 
for most disc stars), and (ii) the stars move outwards much more slowly than their vertical oscillations 
(probably not really the case). We assume that for a given radius, - E^/v - Const., where E. and 
V are the vertical energy and frequency, respectively. From Gauss' law and Poisson's equations, v ~ 
yflnGE, therefore, as stars migrate radially, their vertical frequency changes as v(r) ~ exp(-r/2r^), 
given that the stellar density varies as 2" ~ exp(-r/rrf). To conserve the mean vertical action of a set 
of stars migrating outwards, we then require that their mean vertical energy, {E~) ~ crl, decreases 
as V. Therefore, the vertical velocity dispersion of a stellar sample migrating radially outwards will 
decrease as cr.(r) ~ exp(-r/4rrf). On the other hand, the vertical velocity dispersion in the ambient 
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Fig. 3. Top three panels: initial (dotted blaclc), final (dashed blue), and net (dot-dashed red) vertical velocity 
dispersions, cr-,, cr- f, andzlo"-, respectively, as functions of the initial radius for particles ending up in the green 
bin at the final time. Bottom three panels: the corresponding initial radial distributions of stars ending up in the 
green bin. The bar's CR and OLR are shown by the vertical dotted and solid black lines. Fourth column: time 
evolution of the vertical velocity dispersion profiles for the gSb (top), gSa (middle), and gSO (bottom) discs. 



(non migrating) disc component changes with radius as <Jz{r) ~ yfT{r)h (e.g., eq. 7 by 1261). Assuming 
the scale-height, h, remains constant, we find that cr,(r) ~ exp(-r/2rrf). Thus, while a star's vertical 
energy does decay adiabatically as it migrates outwards, it does not decay as fast as the underlying 
typical vertical energy of disc stars and, therefore, still increases the vertical temperature of the disc 
where it arrives. This is in agreement with our numerical results, where we find that stellar samples 
from the inner disc arrive at 2.5rd with a velocity dispersion slightly larger than the typical values for 
that particular radius (top three panels of Fig. |3]l. 

For the simple model considered, the same argument also holds for the epicyclic energy Ep - kJj,, 
where the radial epicyclic frequency k ~ l/r (valid for a flat rotation curve), so that orbits conserving 
Jp loose energy Ep. This would also suggest that the rates of change in cr^ and cr^ are different for 
the same sample of migrating stars: because in the outer disc at r > r^, v decreases faster than a:, cTr 
is less reduced. This is exactly what we find in plots similar to the upper three panels of Fig. |3] but 
for cr,. (not shown): the outer disc is heated much more radially than vertically due to the migrated 
stars. Consequently, if one were to consider, for instance, an analytic model for radial migration in a 
2D razor-thin disc using Jp as a proxy for extrapolating the full distribution function (including the 
vertical direction), one might be fooled into thinking that a hot thick disc component can be created. 



4 Conclusions 

In this work we have used Tree-SPH N-body simulations to study the effect of internal disc evolution 
on disc thickening. All analyzed models were found to develop strong multiple patterns, resulting in 
very eflicient radial disc mixing. Surprisingly, we found that discs with initial scale-height of ~ 300 pc 
only double their thickness even for bar and spirals much stronger than those of the Milky Way. 
We showed that the reason for this is that stellar samples migrating from the inner disc decrease 
their vertical velocity dispersion. Therefore, although at any given radius the disc heats with time, the 
migrated stars do not contribute significantly to the disc thickness. 

Using a simple galactic disc model, where both the radial, Jp - Ep/K, and vertical, J, = Er/v, 
actions, are conserved, we explained the decrease of velocity dispersion of outward migrators as adia- 
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batic disc cooling. Due to the faster radial decrease of k compared to v, this model predicts that radial 
migration heats discs more radially than vertically, in agreement with the simulations ([16]). While 
stellar samples migrating outwards lose vertical energy, {E-) ~ cP:, their reduced cr^ is still higher than 
what is expected from extrapolating cr^ ~ ^{r). Thus, radial migration should lead to disc flaring. We 
indeed find that scale-height increases with radius in our simulations, however, not significantly so for 
r < 5 scale-lengths. Stellar radial migration also mixes the radial distributions of the actions (7^) and 
the more the disc mixes, the weaker their variation with radius. Therefore, even if the migration 
rate remains constant, we can expect less heating at later times. 

Our simulations did not consider gas accretion. In a more realistic case, the mass increase due 
to continuous thin disc accretion would contract the thickened discs and further decrease the scale- 
heights |2|. Therefore, what we find in this work may be an upper limit on the disc thickening caused 
by radial migration. Considering that the MW thick disc stars are older than 8 Gyr and that thick discs 
are observed at high redshifts, we expect that the 3-4 Gyr of time evolution we follow here is sufficient 
for the problem we address. In other words, if thick discs could be caused by radial migration, then 
that should occur early on, within the first 1/3 of a galaxy lifetime (even much faster if we include 
chemical considerations, e.g., attempt to explain the chemical shifts observed for different abundance 
rations between the thick and thin discs - see [4]). 

For the gSa and gSO models, bulge contribution to the thick disc beyond ~ 2.5 disc scale-lengths 
is only about 1%. For a MW-like bulge (the gSb model) no stars reach this radius (|16|). We con- 
clude that, although radial migration does indeed contribute to the disc thickening, the help from an 
additional mechanism may be needed to account for thick disc formation, such as heating by mergers 
and/or accretion of satellites. At this time the only viable mechanism for purely secular thick disc 
formation remains the heating from turbulent, gas-rich clumpy discs at high redshift ([2 1). 

References 

1. Athanassoula E., 2007, MNRAS, 377, 1569 

2. Bournaud F., Elmegreen B. G. & Martig M., 2009, ApJL, 707, LI 

3. Brunetti M., Chiappini C, & Pfenniger D., 2011, A&A, 534, A75 

4. Chiappini C, 2009, lAU Symposium, 254, 191 

5. Chilingarian I. V., Di Matteo P., Combes F, Melchior A.-L., & Semelin, B., 2010, A&A, 518, A61 

6. Di Matteo, P Combes, F Melchior, A.-L. & Semelin B., 2007, A&A, 468, 61 

7. Famaey, B., Siebert, A. & Minchev, I., 201 1, this volume, ar Xiv;! 1 10-6176 t^l 

8. Foyle K., Rix H.-W. &Zibetti S., 2010, MNRAS 407, 163 

9. Loebman S. R., Roskar R., Debattista V. P, et al. 201 1, ApJ , 737, 8 

10. Masset, F & Tagger, M.,1997, A&A, 322, 442 

11. Meza A., Navarro J. F., Abadi M. G., & Steinmetz M., 2005, MNRAS , 359, 93 

12. Minchev I., Nordhaus J., Quillen A. C, 2007, ApJ, 664, L31 

13. Minchev I., Boily C, Siebert A., et al., 2010, MNRAS, 407, 2122 

14. Minchev I., Famaey B., 2010, ApJ, 722, 1 12 

15. Minchev I., Famaey B., Combes F, Di Matteo P., Mouhcine M., & Wozniak H., 201 1, 527, 147 

16. Minchev I. et al.. In preparation 

17. Patsis P A. & Kaufmann D. E.,1999, A&A, 352, 469 

18. Quillen A. C, Minchev I., Bland-Hawthorn J. &Haywood M., 2009, MNRAS, 397, 1599 

19. Quillen A. C, Dougherty J., Bagley M. B., Minchev I. & Comparetta J., 201 1, MNRAS, 417, 762 

20. Rautiainen P, & Salo H., 1999, A&A, 348, 737 

21. Roskar R., Debattista V. P, Quinn T R., Stinson G. S., & Wadsley J., 2008, ApJ, 684, L79 

22. Scannapieco C, White S. D. M., Springel V., & Tissera P B., 2009, MNRAS, 396, 696 

23. Schonrich R., & Binney J., 2009, MNRAS, 396, 203 

24. Sellwood J. A. & Sparke L. S.,1988, MNRAS, 231, 25P 

25. Shevchenko 1. 1., 201 1, ApJ, 733, 39 

26. van der Kruit, P C. & Freeman, K. C, 201 1, ARA&A, 49, 301 

27. Zaritsky, D. & Rix H.-W., 1997, ApJ, 477, 118 



